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Abstract 

Nonuniform superconducting state due to strong spin magnetism is studied in two- 
dimensional (2D) type-II superconductors near the second order phase transition hue 
between the normal and the superconducting states. The optimum spatial structure 
of the orderparameter is examined in systems with cylindrical symmetric Fermi sur- 
faces. It is found that states with 2D structures have lower free energies than the 
traditional one-dimensional solutions, at low temperatures and high magnetic fields. 
For s-wave pairing, triangular, square, hexagonal states are favored depending on the 
temperature, while square states are favored at low temperatures for d-wave pairing. 
In these states, orderparameters have 2D structures such as square and triangular 
lattices. 

PACS numbers: 74.60.-w,74.70.Kn,74.72.-h,74.76.-w 



Nonuniform superconducting state due to a strong spin magnetic effect, which is called 
Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state, has often been studied, since Fulde & Fer- 
rell and Larkin & Ovchinnikov showed that the upper critical field of a nonuniform super- 
conducting state exceeds the Pauli paramagnetic limit (Chandrasekar-Clogston limit) under 
some ideal conditions.0^ Recently, it is sometimes suggested that organic, high-Tc copper 
oxide, and heavy fermion superconductors can be candidates of the FFLO superconduc- 
tors. This is because they have large upper critical fields so that strong spin magnetism 
is attainable, and also because they can be clean type-II superconductors. 

In particular, the organic and the copper oxide superconductors have large critical fields, 
especially when orbital pair breaking effect is weaken by applying magnetic field in any 
direction parallel to the most conducting plane.§ In addition to this, some features of Fermi 
surface structure arising from the low dimensionality enhances the critical field of the FFLO 
state.iHli Therefore, some compounds in the families of the organic and high-Tc supercon- 
ductors might be good FFLO superconductors. 

However, spatial oscillation of the orderparameter, which characterizes the FFLO state, 
has not been observed yet, while very large upper critical fields were observed in some 
compounds such as (BEDT-TTF)2X and (TMTSF)2X. Those critical fields might exceed the 
Pauli paramagnetic limit In particular, increase of the critical field at low temperatures 
with positive d^Hc/dT^ was observed in (TMTSF)2PF6.0 Such behavior of the critical field 
is very similar to the behavior that is theoretically obtained in low dimensional FFLO 
superconductors. This might suggest the possibility of the FFLO state in this material 
or similar organic compounds, although this is not an evidence of the FFLO state. The 
spatial structure of the orderparameter must be observed in order to prove existence or 
nonexistence of the FFLO state. 

The spatial structure of the orderparameter in the FFLO state is different from that in 
the traditional Abrikosov state, because their origins are quite different. The nonuniformity 
of the orderparameter in the FFLO state is due to the displacement of the Fermi surfaces of 
up and down spin electrons by the Zeeman energy, not due to the orbital motion around flux 
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lines. Gruenberg and Gunther examined the condition for the coexistence of these states 
in a three dimensional system.1 In the two dimensions, the Abrikosov state occurs unless 
the magnetic field is parallel to the conducting layer. However, when the parallel direction 
is approached, its Landau level index increases,0li30 and then the spatial structure of the 
orderparameter approaches to that of the FFLO state.0 We assume a parallel magnetic field 
and concentrate on the FFLO state in this paper. 

In the studies so far, the spatial oscillation of the orderparameter in the FFLO state 
is believed to be in a single direction. Larkin et al. studied an s-wave superconductor 
with spherical symmetric Fermi surfaces, and found that the form A(r) oc cos(q ■ r) is the 
most stable solution among the periodic solutions of the gap equation expanded near the 
second order phase transition field.!! However, other cases such as systems with anisotropic 
Fermi surfaces and those with anisotropic pairing interactions have not been studied yet. 
In particular, two-dimensional (2D) model is significant when one considers the organic 
and copper oxide superconductors. Orderparameters with spatial structures such as the 
square lattice and triangular lattice becomes small in larger area than the orderparameter 
oscillating in a single direction. Thus, the states with such 2D lattice structures would gain 
more spin-polarization energy than the state with the one- dimensional (ID) structures, and 
have chance to occur at high magnetic fields. Hence, such states may occur in the 2D FFLO 
superconductors (including quasi- ID systems), since the critical field of the FFLO state 
remarkably increases at low temperatures. 

In this paper, we study 2D FFLO superconductors with cylindrical symmetric Fermi 
surfaces. We examine s-wave pairing and rf-wave pairing, the latter which is a serious 
candidate of the organic and high-Tc superconductors. Our theory is a straightforward 
extension of the work by Larkin et al. to the 2D systems, finite temperatures, and anisotropic 
superconductivity. We consider the magnetic field parallel to the conducting layers, and 
neglect the orbital pair breaking effect, although weak interlayer interactions are implicitly 
assumed so that the mean field treatment like BCS theory is justified at low temperatures. 

At first, we examine the s-wave pairing. As many authors studied, a gap function of the 
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generalized form 

A(r) = A^e"^-^ (1) 

has the highest second order transition field. The optimum value of |q| is finite for T <T* ^ 
0.56Tio), where Ti°) is the zero field transition temperature. For example, the optimum value 
of |q| is equal to 2h/v-p at T = 0, where h = hqH with the electron magnetic moment hq. In 
the present case, because of the cylindrical symmetry of the system, any linear combinations 

A(r)=^A„e''^-•^ (2) 

m 

have the same critical field, where qm's are any vectors of the optimum magnitude |q^| = q 
parallel to the conducting layer. However, such degeneracy is removed by the nonlinear 
terms in the gap equation, apart from the rotation as a whole: a particular kind of linear 
combinations have a lower free energy than the others, just below the critical field. Here, 
we have taken a finite number of q^'s in the above, because we are considering periodic 
solutions. Near the transition point, the gap equation is expanded as 



log(^)A; = ^[(2 - 6i^)J{eim)AlA^A; 

m 

+(1 - - 5i,^M0im)^l,A_iA*_J, 



(3) 



where 9im is the angle between q; and q^, and 

roo r^TT 



j{eim) = /dc / " -^{iun + ^ + h)-^ 

„ J-oo Jo ZTT 

x(ia;„-^- Vf ■ + h)''^ 



and 



^ j-oo Jo Zn 

X (ia;„ + C + /i + vf • q; + Vf • q^)"^ 
x(iu;„-,^ + /i- Vf -qm)"^- 



(5) 
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We examine periodic solutions 

(a) A(r) = ApF exp(iq • r) 

(b) A(r) = 2Afflo cos(q • r) 

(c) A(r) = 2Asq[cos(gx) + cos{qy)] 

(d) A(r) = Atri[exp(iqi-r) +exp(iq2-r) 

+ exp(iq3-r)] 

(e) A(r) = 2Ahexa[cos(qrr) +cos(q2-r) 

+ cos(q3-r)] 

where qi = 0), q2 = g(— 1/2, \/3/2), and q3 = g(— 1/2, —^/3/2). We refer to the states 
expressed by (a), (b), (c), (d), and (e) as FF state, traditional FFLO state, square state, 
triangular state, and hexagonal state, respectively. 
We define the factor 

|A„|^ = tta^ , (7) 

c 

where a indicates the type of state. For example, a is replaced with FF, FFLO, sq, • • •, for 
the states (a), (b), (c), • • •, respectively. The free energy per unit volume is calculated by the 
formula 

Q - Qo = -JiV(O)^^^ / dh\A{v)\\ (8) 



2 ' ' T, V 



We define the factor ba by 



Q - Qo = -lN(OM\^r. (9) 
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From the gap equation, we have 



1/J(0) 



^FFLO 



2/[J(0) + 2J(7r)] 



'sq 



4/[J(0) + 2J(7r)+4J(7r/2) 



+2J(7r/2)] 



(10) 



3/[J(0)+4J(27r/3)] 



hexa 



6/[J(0) + 2J(7r) +4J(7r/3) 



+4J(27r/3) + 2J(7r/3) + 2J(27r/3)]. 



In the hmit of T ^ 0, we have J(0) = -l/[4/i2(l - ff/^] for g < 1, and J(0) = for g > 1, 
where q = VFq/{2h). Since the optimum value of q is equal to 1 at T = 0, the expansion 
factor J(0) diverges. Thus, the gap equation can not be expanded in the power of A in 
two- dimensions at T = 0. It is easily verified that there is a term proportional to A log A 
in the asymptotic form. 

For finite temperatures, we estimate ba numerically. The results are shown in Fig.|I|. It 
is found that the solution of the form (b), i.e., the traditional FFLO state, is optimum for 
0.24T,(°) <T <T* ^ 0.56T,(°). Below T ^ 0.24T,(°), however, other solutions have lower free 
energies. The triangular, square, hexagonal states are optimum for 0.16Tf ) < T < 0.24T^(o), 
0.05TJ°^ ^ T < 0.16Tj°\ and T < 0.05Tj°\ respectively. This result is plausible because 
when the temperature decreases and the magnetic field increases, node of the orderparameter 
in the real space becomes more favorable owing to gain in the spin-polarization energy. 
Our result is consistent with the result by Buzdin et al., in which the structure of the 
orderparameter is studied near the tricritical point .0 

It is easily verified that the orderparameters of the triangular and hexagonal states have 
spatial structures of the triangular lattice, while that of the square state have a structure of 
the square lattice. The periodicity of the triangular lattice of |A(r)| in the hexagonal state 
is twice that in the triangular state. 
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It is easy to extend the theory to anisotropic superconductivities. The gap function has 
the form 

A(p,r) = 5:A„7(p)e^''"'" (H) 

m 

where 7(p) = a/2(p^ ~Py) ^^"^ d-wave pairing. As Maki et al. studied, the angle between the 
optimum vector q and the x-axis is equal to rni /2 for low temperatures T /T^'^'> < 0.06, while 
7r/4 + n7r/2 for high temperatures 0.06 < T/T^^^ < 0.56, where n is integer.El Since there 
are four equivalent optimum directions of q, we have three candidates to examine, i.e., (a) 
the FF state, (b) the traditional FFLO state, and (c) the square state. 

In this case, a factor [7(p)]^ is introduced into the integrands of eq.(4) and eq.(5). Free 
energies of the above three states are compared in Fig.|^ and Fig.^ Figure ^ shows the 
results for the low temperature phase, which is realized for T /TjP'' < 0.06, while Fig.|^ shows 
the results for the high temperature phase, which is realized for 0.06 < T/Tj:^^ < 0.56. In 
Fig.^], we find that the square state is optimum for T/Tj^^^ < 0.06. In Fig.|^, it is found 
that the square state is optimum for 0.085 < T/T^'^'> < 0.12, while the traditional FFLO 
state is optimum for 0.12 < T/T^") < 0.56. At T/T^ ^ 0.085, the factors and 6sq turn 
negative. Thus, the gap equation expanded up to the third order does not have any solution 
of this type for 0.06 < T/T,(o) < 0.085 below and near the transition temperatures. This 
suggests that a first order transition to the square FFLO state occurs at a field higher than 
the second order transition field, in this temperature region, as discussed below. 

If the gap equation is simply derived from a differentiation of a free energy functional 
expanded in the power of |Asq|, the factor of [Agq]^ in the free energy is negative near 
the second order transition point in this temperature region, since it is proportional to the 
inverse of Ogq and that of 6sq- Thus, the free energy as a function of |Asq| has a minimum 
at nonzero |Asq| at the second order transition point, because there must be higher order 
terms such as a one proportional to |Asq|^, which guarantee that the free energy increases as 
a function of |Asq| where |Asq| is large. This finite value of |Asq| is the optimum solution of 
the gap equation at this point, and this solution must terminate at a higher field. Therefore, 



7 



it is plausible that the actual critical field is higher than the second order transition field, 
and the phase transition is of the first order. 

A similar conjecture as the above holds for three-dimensional systems. As Larkin et al. 
suggested a possibility of the first order transition to the cubic state,i the same possibilities 
for the triangular state and the square state can not be excluded, because the factors for 
these states are negative in the spherical symmetric case. 

Lastly, we breifiy discuss orbital pair breaking effect in three-dimensional FFLO super- 
conductors (including quasi-2D systems). For the ID solution of the traditional FFLO state, 
since the vector q can be oriented in the direction of the magnetic field, the FFLO state can 
coexist with the vortex state if orbital effect is sufficiently weak.! On the other hand, for 
the present 2D solutions, the competition with the vortex states may be a crutial problem, 
since the vortex functions take two out of three digrees of freedom of the coordinate space. 
Lebed argued that if the magnetic field is sufficiently strong, the quasi-2D system can be 
treated as being essentially purely 2D.i In such a case, the competition mentioned above 
does not occur if the magnetic field is accurately parallel to the conducting plane. 

In conclusion, we have examined FFLO superconductors with the cylindrical symmetric 
Fermi surfaces near the second order transition line, and have found that spatial structures of 
the 2D lattices are favored more than the ID structure of the traditional FFLO state, for low 
temperatures and high magnetic fields, both for the s-wave pairing and the (i-wave pairing. 
It is found that the traditional FFLO state, the triangular state, the square state, and the 
hexagonal state appear for the s-wave pairing as the temperature decreases, while the square 
states occur at low temperatures for the (i-wave pairing. If such structures of orderparameters 
are observed experimentally, for example, directly by the scanning tunneling microscope 
(STM) technique, it can be regarded as a strong evidence of existence of a nonuniform 
superconducting state. 
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FIGURES 



Fig. 1. Factors ba in the free energies of the state (a)-(e) for the s-wave pairing. 

Fig. 2. Factors ba in the free energies of the state (a)-(c) for the d-wave pairing, when the an 
between q and x-axis is equal to mr/2. 

Fig. 3. Factors ba in the free energies of the state (a)-(c) for the d-wave pairing, when the an 
between q and x-axis is equal to 7r/4 + mr/2. 
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